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Abstract 

We investigate the realization of a simple solid-state quantum com- 
puter by implementing the Deutsch-Jozsa algorithm in a system of Joseph- 
son charge qubits. Starting from a procedure to carry out the one- 
qubit Deutsch-Jozsa algorithm we show how the iV-qubit version of the 
Bernstein-Vazirani algorithm can be realized. For the implementation of 
the three-qubit Deutsch-Jozsa algorithm we study in detail a setup which 
allows to produce entangled states. 

PACS numbers: 85.25. Cp, 73.23.-b, 03.67.Lx 



1 Introduction 

The growing interest in quantum computation has stimulated an impressive ac- 
tivity in the field of 'quantum hardware'. Numerous proposals to implement 
quantum bits and simple one-bit and two-bit operations have appeared in many 
different areas of contemporary physics research. Yet the possibility to taylor 
a controllable two-state system is by far not enough to do quantum computa- 
tion. From an engineering point of view the design of a quantum computer is 
difficult because of the enormous complexity of requirements, such as the pos- 
sibility to prepare and to measure states easily, a highly flexible setup with a 
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sufficicntly large parameter space that can be controlled, the maintainance of 
coherence etc. Therefore, the touchstone of practical quantum computation is 
the implemcntation of quantum algorithms. 

This goal might appear too ambitious, in particular if one thinks of realizing a 
set of universal gates [El, i.e. a quantum computer which can do any operation. 
However, in order to implement one particular quantum algorithm it is not 
necessary to go as far as this. Since usually one has to deal with a well-defined 
set of input and output states it sufhces to implement just those gates which 
represent the desired operations on applying the gate to the possible input 
states. This renders the task easier than designing gates which represent the 
operations on application to any state. 

Another advantage of the restriction to a particular algorithm is the pos- 
sibility to use more complex operations (for example iV-bit gates instead of a 
sequence of one-bit and two-bit operations) which require less computational 
time. This could be one way to overcome the limits which are set by a small 
decoherence time in a given physical system. 

Surprisingly, there has been comparatively little activity towards the imple- 
mentation of quantum algorithms in real physical systems. To date, quantum al- 
gorithms have been implemented by using liquid-state NMR §, § |, |, § g |, § , 
in atòmic physics jÏ0| and by optical interferometry |ÏÏJ. A solid-state imple- 
mentation of Grover's algorithm has been proposed theoretically |ï^ |. 

The quest for large scale integrability has stimulated an increasing interest 
in superconducting nanocircuits 13, l·|, [Ï5|, [ï(| [ï^] as possible candidates for the 
implementation of a quantum computer. The recent experimental breakthrough 
for Josephson qubits |Ï8|, [Ï9| |2(| is the first important step towards a solid-state 
quantum computer. 

Naturally the question arises whether, at the present technological level, 
it is possible to implement also quantum algorithms in these systems. Here 
we concentrate on charge qubits {T^, |Ï4[ [Ï5| and show how the Deutsch-Jozsa 
algorithm |pT| , ^| and the Bernstein-Vazirani algorithm [p6| can be run 

on a Josephson quantum computer. We analyse the experiment by Nakamura 
et al. Jï^| in terms of quantum interferometry (2Í| and show that it corresponds 
to the implementation of the one-qubit version of Deutsch's algorithm. By 
generalising this idea we show how the iV-qubit Deutsch algorithm, with N < 3, 
can be implemented. 



2 Deutsch-Jozsa algorithm 

Consider the subset of Boolean functions / : {0, 1}^ — > {0, 1} with the property 
that / is either constant or balanced (that is, it has an equal number of 
outputs as ls). The Deutsch-Jozsa algorithm |2Ï], [2^, [23] determines - for a 
given unknown function / - whether the function is constant or balanced. With 
a classical algorithm, this problem would, in the worst case, require 2 N ~ 1 + 1 
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evaluations of / whereas the quantum algorithm solves it with a singlc cvaluation 
by means of the following steps (here we focus on the refined version by Collins 
et al. @, see also Fig. 1). 

(i) All qubits are prepared in the initial state |0), therefore the iV-qubit 
register is in the state |00 . . . 0). 

(ii) Perform an A^-qubit Hadamard transformation Tí 

l*> E \y) > (^ímA (i) 

ye{o,i} N 

where x ■ y = (xx A yx) ® . . . © (xat Ayj\r) is the scalar product modulo two. This 
is equivalent to performing a one-bit Hadamard transformation to each qubit 
individually. 

( iii ) Apply the /-controlled phase shift U / |2^, |24|] 

\x) ^ \x) , (xe{o,in . (2) 



(iv) Perform another Hadamard transformation TL. 

(v) Measure the final state of the register. If the result is |00...0) the 
function / is constant; if, however, the amplitude ot|oo...o) °f the state |00. . .0) 
is zero the function / is balanced. This is because 

«I00...0) =^r E ^ nx) ■ ( 3 ) 
xe{o,i} N 

We note that it is reasonable to identify functions / whose outputs can 
be mapped into each other by a bitwise NOT. For these functions, the gate 
operations in Eq. @ differ merely by a global phase factor (—1). We will use 
the convention /(00 . . .0) = 0. This reduces the number of gates which need to 
be implemented, by a factor of two and does not affect the exponential speed-up. 

In order to implement the algorithm it is necessary to show that each in- 
dividual step (preparation of the state, gate operations, measurement) can be 
realized in a given system. It is well known how to prepare and to measure 
the states in Josephson charge qubits [|Ï3[ |Ï5|, [ïq . Our task is to demonstrate 
that the gate operations corresponding to ah admissible functions / can be per- 
formed with a single device. An important aspect of our proposal is that the 
gate operations are represented in a basis of superpositions of charge states. 

It is interesting to note that according to Ref. |2Íj the Deutsch-Jozsa algo- 
rithm does not involve entanglement for N < 2, i.e. the N = 2 version can be 
realized with uncoupled qubits. On the other hand, the implementation of the 



algorithm for N > 2 involves entanglement |25 and therefore requires a setup 
of coupled qubits. 
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3 Bernstein-Vazirani algorithm 



The Bernstein-Vazirani algorithm [£6], |23| is analogous to the iV-bit Deutsch- 
Jozsa algorithm described in the previous section, with the difference that the 
function / has the form 

f = a-x®b , («,^{0,1}™, í»e{0,l}) . (4) 

The gate in step (iii) is denoted by U a and assumes the form 

\x) ^ (~l)-® b \x) . (5) 

By measuring the register after running the algorithm once one obtains the 
number a in binary representation which classically would require N function 
calls. The Bernstein-Vazirani algorithm does not involve entangled states at 
all ^7j. This can be seen explicitely by rewriting the gate U a as a product of 
single-qubit gates 

N 

Ua = (-l) b H(<rP) a > (6) 

3=1 

where a,j denotes the jth digit of a in binary representation. Here we have 

(i) 

used the Pauli operators a): acting on the computational basis of the jth qubit 
{|0j), and the definition (o>.) :=: 1\ (with the one-qubit identity operator 

Ji). In particular we choose ai' \0j) — +|0j), ai — — 

Apart from the global phase (— l) h the set of gates U a represents the part of 
the Deutsch algorithm with completely separable gates. As the algorithm starts 
with a product state, no entanglement is involved at any step. We note that by 
rewriting the action of the whole algorithm as 

nu a n\oo...o) = (-i) 6 n(^) OJ = C-i) fc IJI«í> 

3 3 

one sees that it leads trivially to the result. In conclusion, it is possible to realize 
the Bernstein-Vazirani algorithm using uncoupled qubits in complete analogy 
with the implemcntation for the one-qubit and two-qubit Deutsch algorithm. 

4 Josephson charge qubits and implementation 
of algorithms without entanglement 

A Josephson charge qubit |ï^, [Ï5| is a Cooper-pair box (see Fig. 3a) which is 
characterized by two energy scales, the charging energy = (2e) 2 /(2C) (here 
C is the total capacitance of the island) and the Josephson energy Ej <C -E c h of 
the tunnel junction. At low temperatures T <C Ej only the two charges states 
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with and 1 excess Cooper pair on the island are important and the system 
behaves as a two-level system with the Hilbert space {|0), |1)} and the one-qubit 
Hamiltonian 

fllq - (^ch/2) {2n x ~l)a z - (E.,/2) a x . (7) 

Here n x = C g V g /(2e) is the offset charge which can be controlled by the gate 
voltage. 

In a recent experiment, Nakamura et al. have measured Rabi oscillations in a 
Josephson charge qubit |Q . In the following we analyse briefly the experiment 
and argue that it can be regarded as the one-qubit implcmcntation of Deutsch's 
algorithm. 

The experiment is shown schematically in Fig. 2. Let us first consider the 
part between the two sudden switchings of the gate voltage. The system has 
been prepared in a superposition of the two states |+), |— ) which are defined as 

l±> = ^= (|0> ± |1)) ■ (8) 

Now a rotation 

exp(i(Ejt/2?í)a z ) (9) 
is performed on the state which results in the final state 

i=( |+) + e~ iE ^ n I-) ) . 

Note that here a z is the Pauli operator with respect to the basis {|+), |— )}. 

After sweeping back the gate voltage, a charge measurement is performed. 
That is, the final superposition is measured in the charge basis {|0),|1)}. In 
particular, for t = 2irh/ ' Ej the outcome of an ideal measurement is |0) while for 
t = nfr/Ej the state |1) is found. 

In a spin-related language the experiment could be described as follows. 
The charge states form the z basis, the charging energy corresponds to the 
magnètic field in z direction, and the Josephson energy corresponds to the field 
in x direction (cf. Eq. (0)). The state is prepared with the magnètic field in z 
direction. Switching the gate voltage suddenly to the degeneracy point flips the 
magnètic field to the x direction. Thus, the spin starts to precess around the 
x axis. After the operation time t the field in z direction is switched on again, 
thus freezing the z component of the spin. The latter corresponds to the island 
charge which then is measured. 

In order to display the analogy between Deutsch's algorithm (see Section 
for one qubit and Nakamura's experiment, we rephrase the algorithm in the 
following way In step (i) and (ii) the symmetric superposition of the states of 
the computational basis {|+), |— )} is prepared. On applying Uf in step (Ui) the 
sign of | — ) in the superposition is changed if the function / is balanced, or it is 
left unchanged. The purpose of the second Hadamard gate is the transformation 
of the superposition to a pure state which is to be measured. 
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Now compare this sequence to Nakamura's experiment. It starts with the 
preparation of the symmetric superposition. The Rabi oscillation corresponds to 
the application of the gate Uf which is implemented by choosing the operation 
time appropriately (see Table I). Instead of performing the second Hadamard 
gate the superposition is measured directly. This is possible because each su- 
perposition corresponds to a charge eigenstate: if / is constant, the charge state 
|0) is obtained while balanced / leads to charge state 

The additional advantage of this procedure is that the second sudden sweep 
brings the Hamiltonian back to a regime where the charge states are approxi- 
mately the eigenstates such that the island charge is frozen. This is important 
since the time of the measurement is much bigger than the intrinsic time scale 
h/Ej of the qubit. 

Thus, Nakamura's experiment with the operation times chosen according to 
Table I can be viewed as the implementation of the one-bit Deutsch algorithm. 
The only difference compared to the sequence of steps in Section II is that 
the single-qubit superpositions are prepared and measured directly, without 



performing additional one-qubit operations 28. . Of course, the essence of the 



algorithm is not affected: instead of rotating the state 'forward' by a Hadamard 
operation we apply a Hadamard rotation 'backwards' to the basis. It is easily 
seen that the coefficient a| 00 ...o) of the charge state |00 . . . 0) still obeys Eq. (||) 
where now the sum has to be taken over all x € {+, —} N - 

The two-qubit Deutsch-Jozsa algorithm can be implemented by using two 
uncoupled qubits The gates a^p implementing 
the balanced functions (the upper index denotes the qubit number) and 1^ ® 

(2) 

I[ for the constant function can be realized in complete analogy to the one- 
qubit algorithm. If, on measuring the qubits, both of them are found in the 
charge state |0) the function / was constant; if at least one qubit is found in 
the state |1), / was balanced. 

Finally, it is also obvious that the Bernstein-Vazirani algorithm can be real- 
ized by applying the procedure described above. According to Eq. (||) one needs 
to implement the operation (ai^) aj for the jth qubit which is straightforward 
by using the entries of Table I (note that the number b in Eq. @ is irrelevant 
for the implementation) . The measurement of the register then yields precisely 
the binary representation of a. 



5 Implementation of algorithms involving en- 
tanglement — the three-qubit Deutsch-Jozsa 
algorithm 

The realization of the three-qubit version of the algorithm is more difficult. 
Apart from the constant function 35 balanced functions need to be implemented. 
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Moreover, the three-qubit algorithm involves gates Uf which produce entangled 
final states. 

The goal is to proceed along the same lines as above, that is, preparation of 
the state |000), sudden sweep of n£ etc. The action of the gates Uf takes place 

in the basis {| + ++), | + H — ), . . . , | )}. That is, the gates operate at the 

degeneracy point rSP = 1/2 of the charge qubits. In order to find efficient ways 
for the implementation we first analyse the functions / and the corresponding 
gates Uf. 

Apart from the constant function and its gate ® ® ij 3 ^ there are 7 
balanced functions for which the gates are separable: ai 1 ) <S> <8 I± , ® 

(2) (3) (1) (2) (3) 

Furthcr thcrc arc 12 balanced functions 
for which the gates factorize into a one-qubit part and a two-qubit part as in 
example 7) below. The other gates entangle all thrcc qubits and can be divided 
into two classes (see example 77) and III))- There are 12 gates of class 77) and 
4 gates of class 777). 

7) \ (l« ® if + a« ® 7Í 2) - J« ® + *W ® af)) ® af 

II) \ (a^) ® l[ 2) ® 7Í 3) - J« ® 7f } ® a^) + ai 1 ) ® a^ ® if + 

+ /í 1) ®ai 2) ®ai 3 )) 

777) i (ai 1 ) ® 7j 2) ® 7j 3) - 7 X (1) ® a| 2 ) ® 7j 3) + 7« ® 7< 2) ® a< 3 ) + 

+ ai 1 )®ai 2 )®ai 3) ) 

All separable single-qubit operations can be carried out in analogy with the one- 
qubit case above. In the following we discuss how the entangling gate operations 
can be achieved. For the realization of these gates a coupling of tunable strength 
between the qubits is required. 

There are various ways to couple charge qubits jïl| [H], |Ï5[ |Ï7). Here we 
investigate coupling via Josephson junctions |2{J. Each qubit island is coupled 
to its nearest neighbour using a symmetric SQUID (see Fig. 3b). 

Assuming that both the jth qubit and the jth coupling junction are tun- 
able by local fluxes $( J ), $^ the Hamiltonian for the iV-qubit system at the 
degeneracy point ns —1/2 reads 

N 

77v q = £ {<(* Ü) ) + ÉgoPoU+V 

-\l/2)4\^)[afa^ + h.c.]} . 
Here JjP is the Josephson energy of the jth coupling SQUID and er± 
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(10) 

= {Ux ± 



ia y )/2. In the limit of small coupling capacitances Cjp <C C^ 7 ') we have 

~ 2(70) ch ' 

(i) 

We will assumc that E£ is negligible. Since in practise the capacitive coupling 
is always present it is necessary to have Jjpfàz: = 0) S> 4E^. Then the 
dynamics of the system approximates the ideal dynamics sufficiently well. 

Considcr now the first and the second qubit coupled by Jjp. By choosing, 

e.g. — E^-p = E^p = ±jj^ = J and the operation time t ~ 0.97(27r/J) we obtain 
an operation similar to a swap gate for the qubits 1 and 2 for which we introduce 
the notation (in the basis {| ++),..., | )}i2) 



/ 











±i 







1 
















1 





V 


±i 












[±12]- n n ï n ■ (") 



By denoting one-bit phase shifts for the jth qubit 

[±i]:=:(j , d2) 

we can write a sequence of operations which gives the two-bit entangling gate 
in example /) above: 

[+1][ -2] — [-12]—*W . (13) 

After suddenly sweeping nP and to the degeneracy, first the one-bit phase 
shifts are performed by keeping J^p = 0. Then J^p is switchcd on suddenly in 
order to do the two-bit rotation. The ai^ rotation can be done at any moment 
since the third qubit is decouplcd from the othcr two. Finally the ni are swept 
back suddenly and the register is measured. Note that the parameters of the 
one-bit and two-bit operations need to be chosen in a compatible way, i.e. the 
local Josephson couplings Ej ' should be the same for the one-bit and two-bit 
operations in order to avoid unnecessary parameter switching. 

There are numerous ways to represent the three-bit entangling gates. At 
least two different two-bit rotations need to be applied. During the second two- 
bit rotation the third qubit has to be 'halted'. This can be done by switching 
off both the Ej and the Jk which couple to this qubit. A possible sequence for 
example II) is 

[+l][-2]-[+13]-[-12] , (14) 

and for example III) 

^V^-I+121-1+231-1+12] . (15) 
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It is interesting to note that the completely entangling gates of class II) 
and can be realized approximately with a single three-qubit operation. In 
Table II and III we list the parameters for the various implementations includ- 
ing estimates for the accuracy of the respective operation. The complete set 
of entangling gates can be obtained from the sequences (13) - ( |Ï5| ) by cyclic 
pcrmutations of qubit numbers (and appropriate sign changes), thereby pay- 
ing attention that the parameter settings are compatible for both one-bit and 
two-bit operations. 



6 Conclusions 

We have presented a possible implementation of the Deutsch-Jozsa algorithm 
in a setup of Josephson charge qubits. While this algorithm for a qubit num- 
ber N > 3 requires entanglement we have demonstrated explicitely that the 
Bcrnstcin-Vazirani algorithm does not involve entanglement for any number of 
qubits. 

Our implementation realizes the algorithms by using only state-of-the-art 
technology. A peculiarity is that the gate operations representing the algorithm 
are carried out in a basis which is different from the one which is measured. This 
helps us to obtain the desired results with a minimum number of operations. 
Thus, one may hope to see the expected behaviour of the system even with the 
relatively small decoherence times which have been achieved up to now . Of 
course, these measured decoherence times refer to a single qubit; at the moment 
it is not clear how much more difhcult it is to observe entangled charge qubits 
experimentally. In fact, the experimental implementation of this proposal may 
serve to study this question in detail. 

Finally, we mention that the methods outlined above can be used to study 
also other interesting problems such as the production and measurement of 
Bell states and GHZ states |3(J. From a practical point of view, it would be 
particularly interesting to find ways to create such states in a 'single shot' with 
one appropriate gate operation. Even though it appears rather difhcult to avoid 
the locality loophole in this kind of setup it is nevertheless a remarkable challenge 
to measure such quantum correlations in a macroscopic system. 
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Table 2: Junction parameters for various realizations of the gates 77) and 777). 
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gate 


implementation 


time t/(2nh/J) 


a |0...0> 

(E$ = 0) 


O|0...0) 

(£«=J/40) 


II) 


sequence (9) 


0.97 (2bit op.) 


2 • 10~ 3 


2 • 10~ 4 


II) 


single operation 


0.80 


7 • 10~ ü 


2 • iï)- 2 


HI) 


sequence (10) 


0.97 (2bit op.) 


3 • 10~ 4 


6 • 10~ 3 


HI) 


single operation 


1.19 


< io- b 


2 • 10~ 3 



Table 3: Gate operation times and accuracy for gate realizations of examples 77) 
and 777). The coeficient a| _ ) is a measure for the fidelity of the operation 
(for an ideal operation a|o... ) = 0, cf. Eq. (||)). The operation time for the 
sequences refers to the time needed for the two-qubit rotations. Single-qubit 
rotations are assumed to be perfect. 
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Figure 1: The sequence of operations to perform the Deutsch algorithm on a 
register of N qubits. According to Ref. |2j| it can be interpreted in terms of 
quantum interferometry. The first Hadamard transformation produces a super- 
position of all possible states. Thus, with the application of the /-controllcd 
gate Uf the outcome of / for all possible arguments is evaluated at the same 
time. The second Hadamard transformation brings all computational paths 
together. 



Figure 2: The qubit is prepared in the ground state |0). After suddenly sweeping 
the gate voltage the system starts Rabi oscillations between the eigenstates of 
the new Hamiltonian |±) = (|0) ± \l))/^/2. After the time t the gate voltage is 
swept back suddenly which freezes the final state; then the qubit is measured. 



Figure 3: a) A charge qubit. The Josephson energy of the junction can be 
controlled by the magnètic flux <£>: Ej(&) — 2£ j cos (ir$>/§o), where £j is the 
Josephson energy of the junctions of the symmetric SQUID and $o = h/ (2e) is 
the flux quantum |Ï5| ]. b) A possible realization of coupled charge qubits. 
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